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ABSTRACT 

We investigate the supercluster shape properties of the all-sky observed 
(Abell/ACO) and simulated (Virgo data) cluster catalogues using an approach 
based on differential geometry. We reliably identify rich superclusters by ap- 
plying percolation statistics to both observed and mock cluster populations, 
the latter being constructed following the observational requirements of the 
Abell/ACO sample extended out to Zmax < 0.114. We apply a set of shape 
diagnostics in order to study in a thorough way the morphological features 
of superclusters with > 8 cluster members. Results demonstrate that filamen- 
tarity is the dominant shape feature of superclusters. On comparing data and 
models, we show that the ACDM (JIa = 1 — = 0.7) model performs better 
than rCDM, which is excluded at a high confidence level, in agreement with 
other recent large scale structure studies. 

Keywords: cosmology: theory - clusters - superclusters: general - large- 
scale structure of universe - Optical: clusters 



INTRODUCTION 



Superclusters of galaxies possess an eminent position in 
the structure hierarchy being the largest mass units we 
observe today. Since they have been seeded by density 
perturbations of the largest scale (~ 100/i~^Mpc), they 
therefore constitute objects with which one can study 
the details of the fluctuations that gave rise to them 
(West 1989; Einasto et al. 1997). Additionally, we can 
obtain fruitful clues regarding the formation and evo- 
lution of the universe on the grandest scales and test 
cosmological models (Bahcall & Soneira 1984; Bahcall 
1988; Frisch et al. 1995). 

Numerous studies have been devoted not only to 
delineate the geometrical pattern of the universe as a 
whole (Zel'dovich, Einasto & Shandarin 1982; de Lap- 
parent, Geller & Huchra 1991) conflrming the picture 
of a well designed web-like network, but also to extract 
statistically complete supercluster catalogues in order 
to analyse their spatial distribution and morphological 
characteristics (Einasto et al. 2001a, b and references 
therein). However, only very recently (Sathyaprakash 
et al. 1998a; Basilakos, Plionis & Rowan- Robinson 2001 
hereafter BPR) was any signiflcance given to cosmo- 
logical inferences from supercluster shape statistics. 
Sathyaprakash et al. (1998a) and BPR have used in- 
frared galaxy samples (1.2 Jy and PSCz respectively), 
whereas we presently consider optical clusters (Abell + 
ACQ). 

A substantial number of geometrical and topologi- 
cal techniques have been applied to observational data 



to explore and assess clustering and superclustering at 
large (> 100/i"^Mpc) distances (Mecke et al. 1994; 
Sahni & Coles 1995; Yess & Shandarin 1996; Dave et 
al. 1997; Kerscher et al. 1997; BPR; Kerscher et al. 
2001a, b and references therein). Nevertheless, an inno- 
vative method was devised recently (Sahni 1998; Sahni 
et al. 1998; Sathyaprakash et al. 1998b) and has been 
used with success, so far, in order to describe in detail 
the global geometrically complex features of large scale 
structure. The flrst time this new scheme was applied 
to astronomical data was in BPR, where it ascertained 
that the prominent feature of the large scale structures 
we see today is filamentarity, as it has also been ob- 
served in N-body simulations of gravitational clustering 
(Sathyaprakash et al. 1998b and references therein). 

The aim of the present work is not to produce a re- 
liable all-sky supercluster list. Such attempts have been 
copiously presented by other authors the last twenty five 
years or so, based either on optical or X-ray (ROSAT) 
cluster data (Bahcall & Soneira 1984; Batutski & Burns 
1985; Postman, Huchra & Geller 1992; Cappi & Mauro- 
gordato 1992; Plionis, Valdarnini & Jing 1992; Zucca et 
al. 1993; Kalinkov & Kuneva 1995; Einasto et al. 2001a, 
b and references therein) . The present analysis will focus 



(i) whether or not superclusters verify the dominance 
of filamentarity as the basic trait of large scale structure 
and 

(ii) if so, whether or not supercluster shape and size 
statistics can be used to test cosmological scenarios. 
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The layout of the paper includes the following 
sections. We first present the observed and simulated 
datasets. In section 3, we discuss useful selection par 
rameters of the observed sample, give a brief but suf- 
ficient report on the method used to investigate super- 
cluster shape properties and comment on its stability. 
In section 4, we derive the morphological parameters 
of the real (Abell/ACO) and mock (Virgo) supercluster 
data and in section 5 we statistically compare data and 
dark matter (DM) models in an attempt to discriminate 
between possible structure formation scenarios. Finally, 
our concluding remarks are deserved in section 6. 



2 THE DATASETS 

2.1 The Abell/ACO cluster catalogue 

In the present analysis we make use of the combined op- 
tical Abell/ACO cluster fist (Abell 1958; AbeU, Corwin 
& Olowin 1989), which contains 2712 and 1364 entries 
in the North and South hemispheres respectively but ex- 
cludes the 1174 supplementary poor southern systems. 
We are utilising a redshift updated version of the cat- 
alogue up to a maximum redshift of amax < 0.114 and 
restrict ourselves to |&| > 30° to avoid severe selection 
biases (cf. sections 2 and 5 of Einasto et al. 2001b) and 
light absorption by the zone of avoidance (ZoA). After 
taking into account the double cluster entries, our fi- 
nal list includes 926 objects (523 AbeU and 403 AGO), 
^ 80% of which (733/926) have measured redsliifts and 
are also in accordance with the above observational lim- 
its. Estimated redshifts are taken from Plionis & Val- 
darnini (1991). 

All the above redshifts are heliocentric and trans- 
formed to the Local Group frame using the latest esti- 
mates for the Sun's galactic coordinates. Redshifts are 
converted to proper distances using a spatially flat back- 
ground cosmology with Ho — lOO/i km s~^ Mpc~^ and 
go = 0.5. Thus, the redshift cutoff corresponds to a lim- 
iting distance J?max < 315h~^ Mpc. 



2.2 The Hubble Volume simulations 

For the purpose of our study it is necessary to use 
the largest up to date cosmological N-body simulations. 
This is the reason for our utilising the Virgo Consortium 
structure formation models (Frenk et al. 2000). Details 
of the simulations have been presented in Colberg et al. 
(2000) and Evrard et al. (2001 and references therein), 
so we only briefly discuss the main points here. 

The Hubble Volume simulations follow the evolu- 
tion of 10® particles in volumes comparable to the whole 
observable universe. We consider two spatially flat struc- 
ture formation models, namely ACDM with Qa = 0.7 
and rCDM (massive r-neutrino contribution). Further 
details of the two models are given in Table 1. The 
lengths of the boxes arc 3/i~^ Gpc and 2h^^ Gpc for the 
two cosmologies respectively. In both cases, the particle 
mass is nip > 2.2 x 10^^ Mq and both models are 
normalised to the abundance of rich clusters at z = 0, 



Table 1. Parameters of the DM models. 



Model 


f2jn 


Qa 


h 


r 


0-8 




rCDM 
ACDM 


1.0 
0.3 


0.0 
0.7 


0.5 
0.7 


0.21 
0.17 


0.6 
0.9 


2;i-iGpc 
3/i-iGpc 



i.e. (78 = 0.55 (Eke, Cole & Frenk 1996) and to 

COBE results. 

Galaxy clusters are specified using a friend-of-friend 

algorithm with linking lengths of order 0.2 and 0.164 for 
tGDM and ACDM respectively. The sizes of simulation 
boxes certify that an appreciable number of independent 
"observers" will be identified within these vast volumes 
and the large number of simulated clusters (> 10*') will 
ensure better statistics when supercluster shape prop- 
erties are measured (see sections 4 and 5). 



3 GEOMETRICAL AND SHAPE 
FORMALISM 

3.1 Selection parameters 

The supercluster catalogues for both observed and 
simulated data are constructed by using a constant 
size neighbourhood radius, i.e. a percolation radius 
(Zel'dovich et al. 1982), a scheme that has been suc- 
cessfully applied to similar kinds of studies. We place 
a sphere of a certain size around each cluster and then 
find all neighbouring spheres having an overlap region. 
In doing so, we join all clusters falling within the spheres 
which have common areas and do this for all clusters in 
the sample. In order to choose the optimal percolation 
radius, we have repeated the procedure by successively 
increasing the size of the sphere. At the end, we identify 
the radius that yields the maximum number of super- 
clusters, which obviously occurs before the percolation 
of the superclusters themselves. The choice of this per- 
colation radius (hereafter pr) has been based on criteria 
similar to those of Einasto et al. (1994; 2001a, b and 
references therein). 

To this end, we plot in Figure 1 the dependence of 
the number of superclusters on pr for a variety of per- 
colation radii for both data and models. Solid line cor- 
responds to the Abell/ACO supercluster data, whereas 
dashed and dotted lines denote ACDM and tCDM re- 
spectively. Note, that this plot accounts for superclus- 
ters with at least 8 members, i.e. with multiplicity 
fc > 8. As the pr increases, the supercluster number 
increases accordingly reaching a maximum at radii be- 
tween 25h~^ Mpc and 32h~^ Mpc. At larger radii, super- 
clusters start connecting with each other and at even 
larger radii, we reach the point of supcrconnectivity 
with only a few (< 5) huge objects pervading the total 
volume available. We, therefore, choose to work with 
Pr = 27h~^ Mpc, which is a value well away from the 
superconnectivity limit (> 33h~^ Mpc). Notice that our 
value is consistent with other recent analyses (Einasto 
et al. 1994; Einasto et al. 1997). 

Another observational parameter we wish to com- 
pute is the number density of the Abell and AGO clus- 
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Figure 1. Number of superclusters N^^i as a function of pr 
for both data and the two DM models. Solid line denotes 
the observational sample, while dashed and dotted lines cor- 
respond to the two models (ACDM, rCDM). Note that the 
superconnectivity limit is similar in all three datasets, oc- 
curring at depths > 33h~^ Mpc. Errorbars (ilc) are shown 
only for one set of data for clarity. 



ter samples as a function of radial depth. This is of 
prime significance since construction of mock cluster 
catalogues will be based on the above parameter. We do 
so by using ten equal volume shells {5V ~ 3 x 10^ h'"^ 
Mpc"^) up to -Rmax < 315/i~^ Mpc and calculate the indi- 
vidual densities of the two observational catalogues sep- 
arately. The corresponding average-over-shells densities 
for the two samples are riAbeii ~ 1-62 x 10~^h^ Mpc~^ 
and nAco ~ 2.55 x 10~^h^ Mpc~^ , giving rise to in- 
tercluster separations of order ~ 39.5h~^ Mpc and 
~ 34/i~^ Mpc respectively. For completeness, we note 
that the number density of the combined catalogue is 
^ 1.82 X 10~^h^ Mpc~^ corresponding to a mean separa- 
tion of ~ 38h~^ Mpc in accordance with other estimates 
of the same population. 



3.2 The shapefinders 

An appreciable number of geometrical and topological 
tools have been developed recently for identifying cor- 
rectly and quantifying the rich texture of large scale 
structure as it appears in large angular or redshift sur- 
veys and N-body simulations (Babul & Starkman 1992; 
Luo & Vishniac 1995; Sahni & Coles 1995; Sahni 1998; 
Sathyaprakash et al. 1998b and references therein). In 
contrast to traditional schemes (low-order statistics) 
which bear no relation to either topology nor shape, 
the shapefinders have been constructed from Minkowski 
functional (volume V, surface area 5", integrated mean 
curvature C and genus G) in an attempt to discrim- 
inate between complex morphologies such as filaments 
(prolate-like objects), pancakes (oblate-like objects), tri- 
axial structures (ribbons) and spheres (see Figure 1 of 
Sahni et al. 1998 and Figure 2 of Sahni 1998 for a 
schematic presentation of the above morphological cat- 
egories) . 

For a particular surface of integration class it 
is well known that we can define the four Minkowski 
functionals which completely characterise the geometry 



Figure 2 . Performance of shapefinders Ki , K2 as a function 
of multiplicity k. Theoretical values (horizontal solid lines) 
are compared with MC computed ones (dashed lines with 
errorbars) for two ideal objects, i.e. a pancake (upper and 
lower lines denoting Ki and K2 respectively) and a triaxial 
object (middle line showing Ki). Errors correspond to the 
±1(T scatter from 100 MC realizations performed for each fc, 
shapefinder and object separately. 



as well as the topology of a compact surface. We, thus, 
introduce the three shapefinders Hi = VS~^, Ti,2 = 
SC~^ and Tis = C having dimensions of length. Based 
on these, we can further define the two dimensionless 
shapefinders Ki and K2 as follows: 



Ki = 



7^2 — Ti-i 



7^2 + Til 
and 



K2 = 



Ti-s — TI2 



(1) 



(2) 



Ti-s + TC2 

where Ki^2 < 1 by definition. 

Therefore, having identified superclusters using the 
specified percolation radius, shape detection is per- 
formed only for those objects having fc > 8 members 
via the moments of inertia method. By fitting the best 
triaxial ellipsoid with parametric form 

r{6, (f)) = Ii sin9 coscfii + I2sm6 sm(l> j + I3COS6 k , (3) 

to real and mock data, we can obtain the three eigenval- 
ues Ji , I2, I3, which directly relate to the three principle 
axes of the fitted ellipsoid. The volume can be then de- 
fined as F = ^ 7i 72 /s and < < 27r, < < TT. 
Subsequently, the latter parameters can further be used 
to derive the shapefinders Hi and Ki. 

The above shape technique characterises the shapes 
of topologically non-trivial cosmic objects according to 
the following classification: 

(i) Pancakes if Ki /K2 > 1 

(ii) Filaments if K1/K2 < 1 

(iii) Triaxial structures if K1/K2 ~ 1 and 

(iv) Spheres if I\ — I2 = I3 and Ki ~ K2 ~ 0. 

Ideal filaments (0,1), pancakes (1,0), triaxial structures 
(1,1) and spheres (0,0) are represented by the four ver- 
tices of the shape plane in the form of the shape vector 
K = {Ki, K2), whose amplitude and direction deter- 
mines the morphology of an arbitrary 3D surface. Fi- 
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Figure 3. Two dimensional whole sky map of the 26 Abell/ACO superclusters (crosses). Open symbols denote the clusters 
associated with superclusters with fe > 8. Note, that we have detected 11 superclusters in the North and 15 in the South. 



nally, for shapes close to spherical, Ki and K2 are very 
small, so that it is meaningful to measure their ratio 

TZ = K1/K2 and evaluate deviations from sphericity 
(for further details see section 4.2 of BPR). 



3.3 Stability of the method 

The present study is based on superclusters with fc > 8, 
as quoted above. This choice has been made on stability 
grounds and is a compromise between adequate statis- 
tics and efficiency of the shape technique (see Plionis, 
Barrow & Frenk 1991; Jaaniste et al. 1998 their sec- 
tion 3.2). Indeed, it has been established that superclus- 
ters shape properties are accurately estimated as long 
as fc > 8 members (Jaaniste et al. 1998 their figure 2). In 
order to investigate the performance of our shape find- 
ing method, we utilise a Monte Carlo (MC) approach. 
In particular, given as input parameters the three axes 
7i, I2 and I3 of an ellipsoid with equation 



222 

X-i Xn I's , , 
— -I -A ^ < 1 

^1 ^2 ^3 



(4) 



for fc > 5, the question that we answer is "what are the 
most probable measured Ki , K2 recovered by our shape 

algorithm?" 

In Figure 2, we compare theoretical predictions 
with MC estimated values of Ki and K2 as a function 
of multiplicity fc. Horizontal solid lines correspond to 
the input values for an ideal pancake (upper and lower) 
and an ideal triaxial structure (middle) and are taken 
from Sahni et al. (1998 their Table 1). Upper (Ki) and 
lower {K2) dashed lines with errorbars denote the de- 
rived figures for the pancake, while the middle one (Ki) 
indicates the triaxial object. For the latter case, we ob- 
serve that our MC computed values always approximate 
extremely well (within 5%) the input ones for all mul- 
tiplicities used. Similar results are also obtained for Ki 
in the case of the pancake (upper lines). Nevertheless, 
it is evident that our K2 estimated figures lie system- 
atically above the input ones, albeit asymptotic conver- 
gence is gradually restored. This systematic trend effec- 
tively means that our algorithm tends to make pancakes 
less planar, a feature which is more pronounced at small 
multiplicities. Exactly the opposite event occurs in the 
case of filaments for the input and MC estimated values 
of Ki, thus producing less elongated filaments than it 
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should. This effect tends to produce a slightly distorted 
shape spectrum of superclusters, keeping, however, the 
shape identity of the sample unchanged (filaments re- 
main filaments and pancakes remain pancakes). 

We have concluded that, for k < 8, our input results 
differ significantly from the MC computed values for 
the study cases of the pancake {K2) and the filament 
(Ki). This is also in line with previous analyses on shape 
determination, thus proving the caveat of the technique 
for superclusters with a few cluster members. 



4 MORPHOLOGICAL PROPERTIES 

4.1 The Abell/ACO supercluster features 

We investigate supercluster characteristics according to 
the observational requirements and the definitions of 
shape diagnostics set in section 3. Taking the optimum 
value of pr = 27/1"^ Mpc for the combined Abell/ACO 
sample, we end up with 26 superclusters with fc > 8. The 
relevant information is presented in Table 2. In Figure 
3, we show a 2D schematic representation of superclus- 
ters (crosses) superimposed onto the related cluster dis- 
tribution (open symbols) for the observed sample. The 
absence of objects in the central stripes of the plot is 
due to the ZoA restrictions that we have placed. Note, 
that only objects with fc > 8 are plotted. 

As is evident from Table 2, there are no spheri- 
cal superclusters in concordance with numerical N-body 
simulations of gravitational clustering and similar stud- 
ies on observed data (Plionis et al. 1992; Sathyaprakash 
et al. 1998a, b; Sahni et al. 1998; BPR). More than haff 
of the systems (54% or 14/26) reveal filamentary struc- 
ture (cf. section 5.2), 5 of which are dubbed as extreme 
filaments (7^ < 0.5). The rest of the objects (12 or 46%) 
are estimated to be pancakes, with 2 of them being ex- 
treme cases {TZ> 2). 

We plot in Figure 4 the results of this study, where 
the ratio TZ as a function of k is shown for all 26 real su- 
perclusters. The equilibrium between filamentarity and 
planarity, however, seems to be perturbed if we extend 
the k infimum from 8 to 10 members. In this case, we 
find 71% (12/17) filaments, while for fc < 10 superclus- 
ters we obtain 64% (7/11) pancakes. We attribute this 
to the ill-defined shapes, as evidenced from the analysis 
of section 3.3 (cf. Figure 2). We thus observe that low 
fc superclusters are mostly compatible with pancakes or 
extreme pancakes (see object at the upper left of Figure 
4), whilst large ones are better described as filaments 
or extreme filaments. Statistically speaking, the larger 
the multiplicity, the more accurately supercluster mor- 
phologies are defined, the more filaments are revealed. 



4.2 The simulated distributions 

The task of constructing mock supercluster catalogues 
based on the two aforementioned DM models (rCDM, 
ACDM) is accomplished by using the selection param- 
eters of the Abell/ACO population presented earlier. 
Since the lengths of the simulated boxes are 2h~^ Gpc 
and 3/i~^ Gpc for the two models respectively, we can 




Figure 4. Ratio 7?. as a function of multiplicity k for the 26 
real superclusters. The vertical dashed line corresponds to 
the fc-cut discussed in text and the horizontal one splits the 
regions of filaments and pancakes. 



easily define several mock cluster distributions as seen 
by a variety of totally independent observers. The num- 
ber of these observers per simulation box can then be 
computed as A'^obs = (ibox/2i?max)^. For rCDM, we 
have defined 27 such observers and 64 for ACDM, each 
of them having exactly the same observational features 
as the initial Abell/ACO population (i.e. number den- 
sity, number of clusters, selection function, geometry 

and i?max). 

We, furthermore, make use of the percolation pa- 
rameter Pi that ma:ximises the mock supercluster num- 
bers for each model (see Figure 1). For rCDM, we have 
computed p^ — 32h~^ Mpc and pr = 28h~^ Mpc for 
ACDM. The latter figures correspond to ~ 21.5 ± 3.5 
superclusters for rCDM and ~ 24±4 for ACDM, where 
means and scatters emanate from the 27 and 64 ob- 
servers respectively. 

It is evident that the choice of Pi both for the mod- 
els and the observational data is essential not only for 
the detection of rich superclusters, but for getting a bet- 
ter handle on cosmic variance effects as well (see next 
section). To this end, we have cross checked our findings 
using the nearest neighbour statistical measure (critical 
radius Ra) for clustered distributions (Peebles 2001) 



Rc 



■ 7 



(5) 



where uis and (n) are the solid angle and mean number 
density and ro, 7 correspond to the correlation length 
and the slope of the correlation function of the sam- 
ple under study. Using a plausible set of parameters 
for the clustering properties (correlation lengths, am- 
plitudes and slopes) of the two models (see section 3.2 
of Colberg et al. 2000) and the average values from the 
literature for the observed sample, we find that our fig- 
ures are very well computed and in excellent agreement 
with those predicted by equation (5) . 



5 SUPERCLUSTER STATISTICS 
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Table 2. List of the Abell/ACO superclusters using = 27h~^ Mpc. The correspondence of the columns is as follows: index 
number, multiplicity, axes of triaxial ellipsoid, distance from us, Ki, K2, their ratio TZ, the morphological classification type 
and finally the supercluster name (see Einasto et al. 2001a their Table Al). F denotes filaments, P is for pancakes and eF and 
eP symbolise extremely elongated or planar superclusters. Note that Ii, I2, I3 and R have units of Mpc. We also report 
that due to our cut in b, only part of the Shapley supercluster has been detected. 



Index 


k 


h 


h 


h 


R 


Ki 


K2 


n 


Type 


Name 


1 


16 


91.6 


26.8 


14.1 


267.9 


0.1405 


0.2622 


0.54 


F 


Pegasus-Pisces + Aquarius 


2 


22 


63.0 


30.6 


25.1 


291.0 


0.0416 


0.0773 


0.54 


F 


Sculptor 


3 


11 


37.1 


20.5 


14.8 


184.1 


0.0454 


0.0644 


0.71 


F 


Pisces - Aries 


4 


9 


31.1 


23.1 


12.1 


129.1 


0.0723 


0.0463 


1.56 


P 


Pisces 


5 


24 


56.7 


23.8 


14.1 


174.1 


0.0945 


0.1429 


0.66 


F 


Pisces - Cetus 


6 


8 


30.8 


20.4 


3.1 


302.4 


0.4551 


0.1194 


3.81 


eP 


Fornax - Eridanus 


7 


18 


84.9 


19.0 


11.7 


117.8 


0.1206 


0.3355 


0.36 


eF 


Sextans + Leo 


8 


8 


34.5 


21.2 


7.7 


204.3 


0.1778 


0.0972 


1.83 


P 


Leo - Sextans 


9 


13 


55.9 


22.3 


12.6 


185.6 


0.1074 


0.1614 


0.67 


F 


Ursa Majoris 


10 


9 


44.8 


22.5 


10.3 


298.0 


0.1356 


0.1218 


1.11 


P 


Leo - Virgo 


11 


19 


63.1 


34.6 


17.9 


225.5 


0.1004 


0.0938 


1.07 


P 


Virgo - Coma 


12 


9 


36.1 


22.1 


7.7 


275.3 


0.1913 


0.1001 


1.91 


P 


Leo A 


13 


16 


42.0 


27.1 


18.2 


295.5 


0.0431 


0.0458 


0.94 


F 


Draco 


14 


1 ^ 

10 


A7 7 


00. u 


J-O.O 


100. ± 


n 1 'i'i'i 

U. J-OOO 


n7i 9 

U.U ( J-Z 


1.0 ( 


p 


Shapley (part) 


15 


16 


55.7 


37.1 


20.2 


184.8 


0.0742 


0.0570 


1.30 


P 


Bootes 


16 


10 


35.0 


14.3 


9.9 


210.6 


0.0713 


0.1336 


0.53 


F 


Corona Borcalis 


17 


15 


56.5 


19.3 


13.2 


116.2 


0.0835 


0.1843 


0.45 


cF 


Hercules 


18 


10 


37.6 


13.7 


10.3 


240.6 


0.0680 


0.1545 


0.44 


eF 


Aquarius B 


19 


9 


27.1 


11.5 


9.6 


168.0 


0.0488 


0.1048 


0.47 


eF 


Aquarius-Cetus 


20 


20 


71.9 


38.7 


28.3 


233.8 


0.0460 


0.0681 


0.68 


F 


Aquarius 


21 


8 


33.1 


18.3 


5.1 


190.5 


0.2757 


0.1317 


2.09 


eP 




22 


46 


79.6 


59.3 


32.4 


189.7 


0.0638 


0.0428 


1.49 


P 


Horologium - Reticulum 


23 


8 


33.9 


19.4 


6.6 


285.2 


0.2069 


0.1141 


1.81 


P 


Caelum 


24 


12 


43.5 


10.9 


8.0 


263.3 


0.0892 


0.2733 


0.33 


eF 


Microscopium 


25 


13 


45.2 


30.7 


12.8 


205.9 


0.1306 


0.0726 


1.80 


P 




26 


8 


31.3 


17.3 


10.3 


277.8 


0.0737 


0.0806 


0.91 


F 


Grus 



5.1 Cosmological implications 

Wc compare the shape and volume (size) spectrum dis- 
tributions between the data and two cosmological mod- 
els via a standard Kolmogorov-Smirnov (KS) statistical 
test. We then compute the corresponding probabilities 
of consistency between models and data (Vks) and place 
these results in Table 3. In order to investigate whether 
the above two statistics can be used as a cosmological 
discriminant, we firstly use the KS test so as to compare 
the two model distributions. It is evident from Table 3 
that the shape spectrum fails to discriminate between 
the two models, probably echoing the Gaussian initial 
conditions conmion to both cosmologies. The volume 
spectrum constitutes a useful tool, since it yields a zero 
probability (~ 10"^") that rCDM and ACDM are being 
selected from the same parent population. 

We display in Figure 5 the results from the shape 
(left plots) and volume (right plots) spectrum. Upper 
panels correspond to the tCDM model and lower to 
ACDM. Hatched areas and open symbols denote the 
observed data and filled circles are for the two DM 
cosmologies. Errorbars in all panels are the ilu val- 
ues based on the independent observers and the vertical 
solid lines in the left mark the transition limits between 
filaments and pancakes. A characteristic worth observ- 
ing is that in both models the percentage of filaments 
exceeds that of pancakes. The latter finding is in accor- 
dance with N-body results (Sathyaprakash et al. 1998a, 



Table 3. KS probabilities ("Pks) of consistency between data 
and models as well as between the two DM models. Note, 
that in this computation, we have used the total number of 
mock superclusters identified within all independent volumes 
for the two models. 



Pair 


Sh. Spectrum 


Vol. Spectrum 


Real - ACDM 


0.682 


0.271 


Real - tCDM 


0.488 


0.001 


ACDM - tCDM 


0.242 


~ 10-19 



b). A clear indication of the discriminative power of the 
volume spectrum test is presented in the right panels 
of Figure 5. It is apparent that the rCDM distribu- 
tion gives the worst fit to the observed data with prob- 
ability of consistency being a mere 0.001. In contrast 
to this, it seems that the volume spectrum of ACDM 
superclusters better represents that of the Abell/ACO 
supercluster sample with corresponding probability of 
consistency being 0.271. This can be explained by re- 
garding the volume spectrum as a natural outcome of 
the different clustering pattern of the two cosmologies 
considered here (see Jenkins et al. 1998; Colberg et al. 
2000 for a comparison between optical data and Virgo 
models) . 

To validate our findings, we further explore the se- 
lection parameters presented in section 3 and test the 
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Figure 5. Statistical comparison between data (hatched ar- 
eas and open circles) and models (filled symbols). Shape 
spectra are plotted in left panels and multiplicity functions 
in right for rCDM (top) and ACDM (bottom) respectively. 



reliability of our methods. On choosing reasonable val- 
ues for pr, the multiplicity infimum and finally a smaller 
set of independent observers, we have recovered that re- 
sults remain remarkably robust. 



5.2 Optical vs Infrared supercluster shapes 

We attempt to compare the shape spectrum of the opti- 
cal (Abell/ACO) and the infrared (PSCz) superclusters 
in order to probe whether or not the prominence of fila- 
mentarity is retained by all-sky data with diverse global 
geometrical and selection properties. For this purpose, 
we have resorted to the recently constructed all-sky 
PSCz supercluster catalogue (BPR) counting ~ 12000 
infrared galaxies. 

We display in Figure 6 the shape distribution of 
infrared and optical data (hatched histograms) for su- 
perclusters with k > 8 (left) and fc > 10 (right). It is ob- 
vious that the filamentarity dominance is confirmed by 
both kinds of data, although it is less prominent in our 
optical sample for k > 8 (cf. Figure 4 and discussion in 
section 4.1). In fact, for k > 10, we observe that infrared 
and optical data yield almost identical values for the 
fractions of filaments (~ 70%) and pancakes (~ 30%). 

A final word of caution is due here. The above argu- 
ment is not to cast any shadow on the results based on 
the use of the nominal k infimum applied to the present 
analysis. Had we used superclusters with k > 10, we 
would have obtained qualitatively the same results but 
run the risk of being too close to the validity limit of 
the KS test. 



6 CONCLUDING REMARKS 

We have explored the shape parameters of superclusters 
identified in observed optical (Abell/ACO) and simu- 
lated (Virgo DM models) cluster populations, up to a 
depth of 315/i~^ Mpc. We have applied a constant size 
percolation radius pr to both data and DM models in 
order to detect superclusters within the volumes of the 
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Figure 6. Shape distributions for Abell/ACO (hatched ar- 
eas) and PSCz superclusters for fc > 8 (left) and fc > 10 
(right). Poissonian errors are only shown for the PSCz sam- 
ple for clarity. 



given cluster samples. Using the selection features of 
the Abell/ACO catalogue, we have constructed several 
mock supercluster samples based on independent ob- 
servers for the two cosmological models (rCDM and 
ACDM). The investigation of supercluster morphologies 
has been based on a differential geometry method de- 
vised by Sahni et al. (1998a). 

Results indicate that filamentary structures dom- 
inate the volumes considered here over pancake mor- 
phologies, as seen in both data and models. This, fur- 
ther supports the idea that filamentarity is the main 
structural feature of large scale structure, since it has 
been reliably detected in all available infrared and op- 
tical samples analysed up to date. No spherical super- 
clusters are detected in accordance with similar obser- 
vational analyses and theoretical predictions from large 
N-body experiments (Sathyaprakash et al. 1998a; BPR). 

A standard KS statistical test of consistency be- 
tween data and the two models on the size distributions 
has clearly given a preference to the ACDM {Qa ~ 0.7) 
cosmology. On the contrary, tCDM is rejected at a high 
(99.9%) significance level. Finally, the shape spectrum is 
insensitive to the two cosmologies, reflecting their com- 
mon Gaussian initial conditions. The latter findings are 
in concordance with a similar analysis based on infrared 
PSCz supercluster properties. 
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